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A model of hard-core bosons and spin-1 sites with single-ion anisotropy is proposed to approximately 
describe hole pairs moving in a background of singlets and triplets with the aim of exploring the 
relationship between superconductivity and antiferromagnetism. The properties of this model at 
zero temperature were investigated using quantum Monte Carlo techniques. The most important 
feature found is the suppression of superconductivity, as long range coherence of preformed pairs, 
due to the presence of both antiferromagnetism and S z = ±1 excitations. Indications of charge 
ordered and other phases are also discussed. 
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The mechanism of pairing and the establishment of 
long-range superconducting (SC) coherence are still ceru 
tral issues in the theory of high-T c superconductivity.EJ 
Although strong electronic correlations are the essential 
component of most proposed scenarios for the SC phase 
and its nearby antiferromagnetic (AF) Mott insulating 
phase, important experimental results such as the res- 
onant peaka have not been satisfactorily explained by 
such scenarios. Part of the theoretical limitations are 
related to the enormous difficulty in studying a micro- 
scopic model like the t-J model, which is the simplest 
model that describes the dynamics of holes in an AF 
background. From the numerical point of view, the 
main difficulty consists in reaching large enough clusters. 
This problem is really critical in the presence of inho- 
mogeneities, like the well-known stripeslll and the ones 
that have more recently become the center of intensive 
research^ 

In this article, we propose and analyze a highly sim- 
plified effective Hamiltonian in two dimensions (2D) to 
study the interplay between superconductivity and anti- 
ferromagnetism. Our goal is to describe the movement 
of boson hole pairs in a sea of magnetic excitations. This 
model implies a reduction in the Hilbert space with re- 
spect to the t-J model, which is convenient for exact 
diagonalization calculations, and the elimination of the 
"minus sign" problcmcl which severely inhibits the ap- 
plication of quantum Monte Carlo (QMC) techniques to 
fermions in 2D. Such a model could be useful to study 
experimental features like the ones mentioned earlier and 
inhomogeneous states in which there are phase separated 
SC and AF regions. 

Although I shall not attempt a derivation of the present 
model from a more microscopic one like the t-J model, I 
shall present some heuristic arguments to guide the phys- 
ical interpretation of the results obtained and shown be- 
low. Let us start with the t-J model on the square lattice. 
A coarse-grained Hamiltonian can be obtained by map- 
ping two nearest neighbor (NN) sites (a "dimer" ) jof the 
original lattice onto a site of the effective lattice.El Due 



to the constraint of no double-occupancy in the starting 
model, there are nine states in each of the coarse-grained 
sites. Since our purpose is to study the interplay between 
pairs and magnetic excitations, the states corresponding 
to only one hole in the original dimers are eliminated. 
Hence, we are left with five states per site which cor- 
respond to a singlet (S = 0), a triplet (S = 1, and 
S z = —1,0,-1-1), and a hole pair. Finally, the singlets 
play the role of the vacuum, and the hole pairs are de- 
scribed by bosons .El All these states are by construction 
hard-core entities. Of course, this coarse-graining proce- 
dure here outlined has been employed many times, par- 
ticularly in the context of resonant valence bond (RVB) 
theories. By restating the interactions between dimers, 
it has been shownH that the long-range AF order is re- 
covered in the undoped case. 

The physical underlying scenario is then one in which 
hole pairs move in a "soup" of singlets and triplets. The 
same states pet, site can be also found in "projected" 
SO (5) modelsLrli!] but in these models the interactions be- 
tween them are essentially dictated by symmetry require- 
ments rather than by microscopical considerations (see 
Fig. [I] below). An extreme case with no triplet excita- 
tions, i.e., a model of pairs as hard-core bosons, has .been 
extensively considered to describe superfluid phases .IlirMl 

However, the model studied in the present work con- 
tains a further simplification and hence it should be con- 
sidered as a "toy model" of that scenario of pairs moving 
in a singlet-triplet soup. This simplification consists in 
considering doublets instead of triplets. This is achieved 
by assuming that all sites are occupied by triplets and 
pairs, and by adding to the Hamiltonian a single-ion 
anisotropy term (third term in Eq. (|I]) below) with a cou- 
pling constant A. For A ^ 0, the triplet states split into 
S z = ±1 doublets and S z — singlets. Of course, the 
total spin is no longer a good quantum number. Let us 
introduce D = J2i(^i) 2 > the number of doublets. For 
A = there are only triplets in the model and then 
D goes to 2/3 of the number of triplets. For A — > oo 
there are only singlets and then D — » 0. Assuming that 
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this quantity D evolves continuously between both limits 
then one could consider it as the analog, in this simplified 
model, of the number of triplet excitations in the singlet- 
triplet soup earlier discussed. In this sense, the single-ion 
anisotropy A would correspond to the chemical potential 
in the related SO(5) models. 

The exchange interactions included in our model are 
those resulting from a single exchange in the original 
spin-1/2 t-J model, as shown in Fig. [I]. The exchange 
interaction shown in Fig. |](b) is not included in the pro- 
jected SO (5) Hamiltonian. On the other hand, in that 
Hamiltonian, there is a term corresponding to the one 
shown in Fig. |l|(c). However, this interaction has no mi- 
croscopic origin in a t-J like model because it does not 
conserve the total S z and hence it is not included in our 
Hamiltonian. Then, the effective Hamiltonian here pro- 
posed is£j 
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where Cj :S , cj s are annihilation and creation operators of 
s = S z = 0, ±1 spins, bi, b\ are annihilation and creation 
operators of hole pairs, = b\bi. The exchange term 
is just the spin-1 Heisenberg term, and it captures the 
corresponding interactions shown in Fig. ^ in a simple 
way. 
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FIG. 1. Microscopic origin of exchange interactions. 
Dimers are indicated with bold lines. The exchange in- 
teraction between electrons in different dimers are shown 
with dashed lines. In terms of the total S z of each dimer, 
case (a) corresponds to (0, +1) «-> (+1,0), and (b) to 
(0,0) <-> (-1,+1). The change (c), (0,0) <-» (+1,+1) has 
no microscopical origin. 

The hopping term between pairs and doublets (s = ±1 
in the first term of Eq. (|l|)) induces the minus sign prob- 
lem in the quantum Monte Carlo simulations and hence 
we shall not include it in most of the following study. 
The analogous term of hopping between pairs and triplets 
are also not included in the projected SO (5) Hamilto- 
nian. Preliminary results obtained including this term 
show essentially the same qualitative features shown be- 
low for the parameters considered. Although in principle 
A should be determined by the internal dynamics of the 



original t-J model, in our effective model this is a pa- 
rameter which we vary freely, in the same way as the 
spin chemical potential in the projected SO (5) model. 
Finally, our model contains a nearest neighbor Coulomb 
repulsion between pairs,, to prevent phase separation in 
the low density region .Ej 

We adopt t as the unit of energy. From previous stud- 
ies of the hard-core boson model we adopt V = 3. It is 
not simple to determine a priori the ratio J/t. J should 
be of the order of the exchange coupling constant of the 
original Hamiltonian, although in some cases it may in- 
volve longer than NN interactions, t ~ tot' Q /Ai> where 
to (t' ) is the NN (longer than NN) hopping amplitude 
in the original t-J model, and Af, the cost of breaking a 
hole pair in the intermediate state. This ratio could be 
determined eventually a posteriori by e.g. matching the 
energies of small clusters of the original and the effective 
model. In the present work, we adopted a reasonably 
small value J/t — 0.3. We expect that the effect of a 
different value would just amount to quantitative shift- 
ing the phase boundaries without introducing qualitative 
new features. Of course, a more exhaustive parameter 
study should be done in the future. 




FIG. 2. Relative number of sites with S z — ±1 for various 
values of the single-ion anisotropy A. The stars correspond to 
the full Hamiltonian Eq. (|l|) at T = 0.2 (in this case a typical 
error bar is shown). 

All results shown below, except otherwise stated, cor- 
respond to the model Eq. ([!]) excluding s — ±1 in the 
hopping term, and they were obtained by QMC tech- 
niques (conventional world-line algorithm) on the 8x8 
cluster with periodic boundary conditions. The temper- 
ature T was varied between 0.2t and O.lt, i.e., consid- 
erably below the Kosterlitz-Thouless temperatures sepa- 
rating normal and SC phases found in previous studies 
on related modelsJlj'E3 The Trotter number was kept at 
a standard value, At = 0.1. The results were finally ex- 
trapolated to zero temperature with an exponential law. 
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The error bars are about or smaller than the size of the 
symbols used, except otherwise stated. 

The dependence of D/N (N: number of sites of the 
cluster), i.e. the relative number of S z = ±1 sites is 
shown in Fig. |as a function of the pair density x and 
for various values of A > 1. For A = 8, D is negligible 
and hence this value of A is essentially the A — > oo limit, 
i.e., harcLcpte pairs moving in a vacuum played by the 
singlets. lilfEj The continuity of D as A is increased from 
zero, at a fixed x is not a trivial problem. It is obvious 
that for J = there are no terms in the Hamiltonian 
to compensate the cost of the single-ion anisotropy term 
and hence there is a discontinuity in D as soon as A takes 
a nonzero value. The possibility of a finite critical value 
of J, J cr (x), below which D is discontinuous with A will 
not be examined in the present work. It should also be 
noticed that, after dropping the s = ±1 contributions 
to the first term of Eq. ([!]) , another source of anisotropy 
appears between the S z = 0, ±1 components. This is 
reflected in D as it can be seen in Fig. ^ for A = 2. 

In order to determine partially the phase diagram of 
this model in the A — x space, we compute the follow- 
ing quantities. In the first place, charge and spin cor- 
relations, C(r) = (n n r ), S(r) = (S§S"), respectively, 
together with their Fourier transforms C(k) and S(k), 
i.e., the charge and spin static structure factors. In ad- 
dition, we compute also the staggered correlations, i.e., 
including a factor (—l) x+y in the sums of the previous 
expressions. 




FIG. 3. Staggered spin correlations at the maximum dis- 
tance (open symbols) as a function of pair density and for 
various values of A. The charge structure factor at (tt, tt) is 
shown for A = 8 (circles) and A = 1 (stars). In the inset, 
the charge correlations vs. distance are shown for A = 1 at 
x = 0.03125 (multiplied by 10) and x = 0.25. 

In Fig. ||, the staggered spin correlations at the maxi- 
mum distance r, are shown as a function of pair density 
and for several values of A. They indicate strong AF 



correlations as A — > 0, particularly as x — > 0. i.e. as the 
pure spin-1 Heisenberg model is approached. The deter- 
mination of the long-range character and of the boundary 
of the AF region as a function of A would imply a finite 
size scaling, which is out of the scope of the present study. 
The charge structure factor has a peak at (tt, tt) for all 
the parameter space (x, A) examined except at low den- 
sities. As shown in the Fig. ||, this structure factor is 
very weakly dependent with A, and hence the so-called 
"checkerboard solid" phase found in the hard-core bo- 
son limitliil might extend down to A = at x ~ (see 
discussion below). 

For pair densities lower than k 0.15 and for A 
we have detected indications of incommensurate charge 
ordering with qic = (tt — S, tt — 6) below x — 0.06, and 
qic = (tt — 5, tt) (and symmetry related points) above 
it, with 5 varying with pair density. In these cases, the 
charge structure factor is more than two orders of mag- 
nitude smaller than for (tt, tt) at x = 0.5. 

In Fig. ||(a) the total energy per site vs. pair density 
for various A is shown. The region x > 0.4 apparently 
presents a negative curvature which corresponds to phase 
separation (PS), although probably larger clusters would 
be needed to confirm this result. In this region, PS was 
found for A — > oo, 0-39 < x < 0.5 in the bulk limit using 
alternative criteria.E^ 




FIG. 4. (a) Energy per site, and (b) kinetic energy per site 
as a function of pair density and for several values of A. 

Consistently with this indication of PS, pair-pair cor- 
relations (inset of Fig. ||) show a tendency of the holes 
to be at short distances as the density increases. This 
tendency is superimposed to a strong alternation corre- 
sponding to the peak in the charge structure factor at 
(tt, tt). The second most important peaks correspond to 
QIC — {tt — S, tt), S — 7r/4 (and symmetry related points), 
and there are no indications of pair clustering in a com- 
pact region (which would correspond to q = (5,6)). 

To determine the presence of a superconducting phase, 
considering the hard core bosons as tightly-bound Cooper 
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pairs, we compute the equal-time current-current corre- 
lations: 

tfaa(r) = -j^&Wi«(< + *)> (2) 

i 

where the current operator along direction a (a = x, y) 
is: 

j a (i) = it^2(bl +a Ci+a,s4,s b i - h - c -) ( 3 ) 

s 

We adopt as the SC order parameter the correlation be- 
tween the total current crossing a border and the current 
on a reference bond at the maximum distance: 

xsc = Y. c vy^ L l 2 ) ( 4 ) 

X 

The signs of Cyy (r) (y is the vertical direction) are shown 
in Fig. ||. At the maximum xsc (Fig- ||( a ))j tne pattern 
shows the presence of current lines extended over all the 
cluster, while for xsc ~ 0, x ~ 0.5 (Fig. j|(b)), the cur- 
rent loops are of the size of a single plaquette. 




FIG. 5. Signs of C yy (r) (y is the vertical direction), (a) at 
the maximum xsc and (b) for xsc ~ 0, x ~ 0.5. The arrows 
indicate the reference bond 

The most important results obtained with the present 
model concern the superconducting order parameter 
shown in Fig. ^. xsc and Cyy{v) at the maximum dis- 
tance along y have a similar behavior. Even though there 
are important error bars, there are two apparent qualita- 
tive features. In the first place as A is reduced (A = 8 is 
representative of the A — > oo behavior) the onset of the 
SC phase is pushed to higher values of pair density. This 
feature may be a consequence of AF order at larger den- 
sities although, as said before, a finite size extrapolation 
should be done to determine the phase boundaries. Sec- 
ond, the intensity of xsc is reduced as A decreases, i.e. as 
more S z = ±1 sites are present in the system. In other 
words, there is a maximal SC signal in the hard-core 
boson model and this signal is reduced as the isotropic 
spin-1 model is approached. 

It is also instructive to compare the behavior of xsc 
with the one of the kinetic energy per site shown in 
Fig. ||(b). Both xsc and (K)/N decrease with A for a 
given x, although this effect is more intense in the former 
quantity. The suppression of both quantities near 




FIG. 6. Superconducting order parameter xsc as a func- 
tion of pair density and for several values of A. The cur- 
rent-current correlations C yy at the maximum distance along 
y are shown in the inset as a function of pair density and for 
A = 1 and 8. Typical error bars are shown. 

x = 0.5 may be due to the charge localization in the 
(7r,7r) solid phase. Actually, since xsc and C(n, it) are 
sizable for x > 0.4, this suggests a coexistence of SC and 
the checkerboard solid in this .region, which corresponds 
to the PS region found aboveJlj On the other hand, the 
suppression of SC at low densities has no counterpart in 
the behavior of the kinetic energy. The characteristics 
of this region are hence indicating the presence of a AF 
Mott insulating phase. 

0.3 | ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 




FIG. 7. Superfluid density as a function of pair density and 
for several values of A. 

Finally in Fig. [t] we show the superfluid density which 
is a quantity more conventionally studied in the coa-tcxt 
of hard-core bosons. This quantity is computed asta 
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As = / (ja{T)j a )dT 
JO 



(5) 



where r is the imaginary time, (3 — 1/T and j a {T~) = 
J2ija(h T )- As it can be seen in this Figure, p s has the 
same qualitative behavior as the kinetic energy per site. 
Hence, as discussed above, it differs with xsc m the be- 
havior at low pair density. 

In summary, a model of pairs as hard-core bosons in- 
teracting with spin-1 sites has been proposed to study 
the interplay between superconductivity and antiferro- 
magnetism. Using quantum Monte Carlo techniques, its 
main features as a function of pair density and single- 
ion anisotropy have been determined. Even in this "toy" 
model some features, which could have some relevance 
to high-T c superconductivity, have been observed. The 
most important feature is the suppression of supercon- 
ductivity, as long range coherence of preformed pairs, due 
to the presence of both antiferromagnetism and S z = ±1 
excitations, which in this model represent triplet excita- 
tions in more realistic models. Indications of incommen- 
surate charge ordering and phase separation have also 
been observed. This model certainly deserves further in- 
vestigation. In particular, the dynamical correlations be- 
tween doublets and pairs should be examined. Finally, 
the proposed model is relatively simple and hence it can 
be expected to describe other physical systems. In fact, 
the model can describe a variety of spin-1 systems with 
a layered structure like L^-zSr^MiC^ (Ref. |l6|) or with 
two-dimensionally coupled chains £3 Hole pairs could ap- 
pear in the former compound if there is an on-site attrac- 
tion present. 
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